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ABSTRACT 

We study the clustering properties of X-ray sources detected in the wide area 
(~ 2deg^) bright, contiguous XMM-Newton/2dF survey. We detect 432 objects 
to a flux limit of 5 x 10^"'^^ ergcm~^s^^ in the soft 0.5-2 keV band. Performing 
the standard angular correlation function analysis, a ~ Scr correlation signal 
between and 150 arcsec is detected: w{9 < 150") ~ 0.114 ± 0.037. If the 
angular correlation function is modeled as a power law, w{0) = {0o/O)'^~^ , 
then for its nominal slope of 7 = 1.8 we estimate, after correcting for the 
integral constraint and the amplification bias, that 60 — 10.4 ± 1.9 arcsec. 
Very similar results are obtained for the 462 sources detected in the total 
0.5-8 keV band {9o ^ 10.8 ± 1.9 arcsec). 

Using a clustering evolution model which is constant in comoving coor- 
dinates (e = —1.2), a luminosity dependent density evolution model for the 
X-ray luminosity function and the concordance cosmological model (fim = 
1 — flA — 0.3) we obtain, by inverting Limber's integral equation, a spatial 
correlation length of Tq ^ 16 h^^ Mpc. This value is larger than that of previ- 
ous ROSAT surveys as well as of the optical two-degree quasar redshift survey. 
Only in models where the clustering remains constant in physical coordinates 
(e = —3), do we obtain an Tq value (~ 7.5 Mpc) which is consistent with 
the above surveys. 

Finally, comparing the measured angular correlation function with the 
predictions of the concordance cosmological model, we find for two different 
bias evolution models that the soft X-ray sources at the present time should be 
biased with respect to the underline matter fluctuation field with bias values 
in the range (which depends on the biasing model used): 1.9< 2.7 for 
e = -1.2 or 1< 6o< 1-6 for e = -3. 



Keywords: galaxies: clusters: general - cosmology: theory - large-scale 
structure of universe 



1 INTRODUCTION 

Active Galactic Nuclei (AGN) can be detected out to 
high redshifts and therefore, study of their clustering 
properties can provide information on both the large 
scale structure of the underlying matter distribution and 
its evolution with redshift. At optical wavelengths the 
2dF QSO redshift survey (2QZ; Groom et al. 2000) com- 
prising over 25 000 optically selected QSOs in the range 
z ~ 0.3 — 3 has provided tight constraints on the spa- 
tial distribution of powerful AGNs (Groom et al. 2001; 
Groom et al. 2002). A striking result from this survey 
was that the clustering properties of QSOs are compa- 
rable to those of local galaxies. Moreover, when studied 
as a function of redshift the clustering of these sources 
was found to be constant out to z ~ 3. 



Optically selected AGN catalogues however, are be- 
lieved to miss large numbers of dusty systems and there- 
fore, provide a biased census of the AGN phenomenon. 
X-ray surveys, are least affected by dust providing an 
efficient tool for compiling uncensored AGN samples 
over a wide redshift range. From the cosmological point 
of view an interesting question that remains to be ad- 
dressed is how the X-ray selected AGNs trace the under- 
lying mass distribution and whether there are any dif- 
ferences with optically selected samples. Despite the im- 
portance of X-ray selected AGNs, their clustering prop- 
erties remain poorly constrained. Early studies with the 
Einstein and the ROSAT satellites have produced con- 
tradictory results. Boyle & Mo (1993) used low red- 
shift AGNs detected in the Einstein Medium Sensitivity 
Survey (EMSS; reference) and found only a marginally 
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significant clustering signal at scales < 10 ft ^ Mpc. 
Vikhlinin & Forman (1995) combined archival ROSAT 
observations totaling 40 deg^ and detected, for the first 
time, a statistically significant clustering signal using 
angular correlation function analysis. Their results sug- 
gest a clustering length consistent with that of opti- 
cally selected QSOs. Akylas, Georgantopoulos & Plionis 
(2000) used the ROSAT All Sky Survey Bright Source 
Catalogue to explore the clustering of nearby AGNs. 
They estimate Fo = 6.5±1.0 h^^ Mpc, which is also sim- 
ilar to nearby galaxies and the 2QZ survey results. Con- 
trary to the studies above that are based on an angular 
correlation analysis, Carrcra ct al. (1998) used rcdshift 
information to measure the spatial correlation function 
of X-ray sources in the ROSAT Y)ecY> (Georgantopoulos 
et al. 1996) and RIXOS (Mason et al. 2000) surveys. 
They detect only a marginally significant clustering sig- 
nal and argue that their results suggest that the X-ray 
population is more weakly clustered than optically se- 
lected galaxies or AGNs. Recently, MuUis et al. (2004) 
using the ROSAT North Ecliptic Pole (NEP) survey of 
relatively local X-ray selected AGNs, found a spatial 
correlation length of To ~ 7.4 ± 1.8 Mpc within the 
concordance cosmological model. 

The new generation Chandra and XMM-Newton 
telescopes have extended the studies above to the hard 
(2-8 keV) spectral band. Yang et al. (2003) used Chan- 
dra observations and argued that hard (2-8 keV) X-ray 
selected sources have large variance (strong clustering) 
and are most likely associated with high density regions. 
These authors also find that X-ray sources selected 
in the soft (0.5-2 keV) energy band are less clustered 
(about 1 dex) than hard ones. Recently, Basilakos et al. 
(2004) applied an angular correlation function analysis 
to hard X-ray selected sources detected in the wide area, 
shallow XMM-A'^ewton/2dF survey. They find a strong 
signal consistent with a spatial clustering length in the 
range ro ~ 10 — 19 Mpc (in the concordance cos- 
mological model). This also suggests that hard X-ray 
sources could trace the high density peaks of the under- 
lying mass distribution. 

In this paper we further explore the clustering prop- 
erties of the X-ray population exploiting the high sen- 
sitivity and the large field-of-view of the XMM- iVewton 
observatory. In particular, we extend the Basilakos et 
al. (2004) clustering study to sources detected in the 
soft (0.5-2 keV) and the total (0.5-8 keV) spectral bands 
of a wide area (as 2deg^), contiguous 'KMM-Newton 
survey (XMM-A^eu)ton/2dF survey). Our study provides 
the first constraints on the clustering properties of the 
sources in the above spectral bands using the XMM- 
Newton. Furthermore, we model our X-ray source clus- 
tering and its evolution in an attempt to derive their 
present time bias with respect to the underline mass 
fluctuation field. 

The structure of the paper is as follows. The X- 
ray sample is presented in Section 2 and the angular 
correlation function analysis is discussed in Section 3, 
while the spatial clustering predictions arc presented 
in section 4. Section 5 outlines the models used to in- 
terpret the angular correlation function results and the 
theoretical interpretation of the X-ray source cluster- 
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Figure 1. Left panel: The soft band (0.5-2 kcV) differen- 
tial number counts from the present survey in comparison 
with the best fit double power-law to the counts from Baldi 
et al. (2002; dashed line). Right panel: The total band (0.5- 
8 keV) differential counts from the present survey in compar- 
ison with the results from Manners et al. (2003; dashed line) 
spanning the flux range 10~^^ — 8 X 10~^* ergs"'^ cm"^. At 
bright fluxes extrapolation of the Manners et al. (2003) best 
fit relation overestimates the observed X-ray source surface 
density. The continuous line is our best fit relation to the 
0.5-8 keV counts from the XMM-7Ve«)4on/2dF survey using 
the double power-law described in the text. 

ing. Finally, we draw our conclusions in section 6. Here- 
after and wherever necessary we will assume the con- 
cordance cosmological model (unless stated otherwise), 
ie., ^A = 1, = 0.3, Ho = 100 hkm s"^ 

Mpc"^ (Spergel et al. 2003; Tegmark et al. 2004) with 
h ~ 0.7 (Preedman et al. 2001; Peebles and Ratra 2002 
and references therein) and baryonic density parame- 
ter n^h^ ~ 0.02 (cf. Olive, Steigman & Walker 2000; 
Kirkman et al 2003). 



2 THE SAMPLE 

The X-ray data used in this study are from the XMM- 
Newton/2dF survey. This is a shallow (2-10 ksec per 
pointing) survey carried out by the XMM-Newton near 
the North Galactic Pole [NGP; RA(J2000)=13''41'"; 
Dec.(J2000)=00°00'] and the South Galactic Pole [SGP; 
RA(J2000)=00''57'", Dec.(J2000)=-28°00'] regions. A 
total of 18 XMM-Newton pointings were observed 
equally split between the NGP and the SGP areas. A 
number of pointings were discarded due to elevated par- 
ticle background at the time of the observation result- 
ing in a total of 13 usable XMM-Newton pointings. A 
full description of the data reduction, source detection 
and flux estimation are presented by Georgakakis et al. 
(2003, 2004). 

Here we use the soft (0.5-2 keV) and the total 
(0.5-8 keV) band catalogues of the XMM-Newton/2dF 
survey. We only consider sources at off-axis an- 
gles < 13.5 arcmin. The two samples comprise 432 
and 462 sources respectively above the 5a detection 
threshold. The limiting fluxes are fx (0.5 — 2) = 
2.7 X 10-^^erg s-^cm-2 and /x(0.5 - 8) = 6.0 x 
10~^^ erg s~^cm~^. The sensitivity of the XMM-Newton 
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degrades from the center to the edge of the field of view 
(vignetting) and therefore the limiting flux varies across 
the surveyed area. We account for this effect by con- 
structing sensitivity maps giving the area of the sur- 
vey accessible to point sources above a given flux limit. 
In the 0.5-8 keV band about 10 per cent of the sur- 
veyed area is covered at the flux /x(0.5 — 8keV) = 
10"^'' ergs~^ cm~^. This fraction increases to about 50 
per cent at fx{0.5 - SkeV) = 2 x lO"^"* erg s"^ cm"^. 
In the soft band about 10 and 50 per cent of the 
total area is covered at the flux fx (0.5 — 2keV) = 
3.5 X 10"^^ ergs"^ cm"^ and 5 x 10"^^ ergs"^ cm"^ re- 
spectively. 

Unfortunately, the identification of our sources in 
the 0.5-2 keV band remains largely unknown since op- 
tical spectroscopy is not available for the large majority 
of them. However, from other surveys in the same band 
and of similar depth, we know that the vast majority of 
sources are associated with AGN. For example, among 
the 50 soft X-ray selected sources in the ROSAT Lock- 
man Deep Field (Schmidt et al. 1998), reaching a flux 



depth of fo.5- 
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erg cm ^ s ^,65 and 15 per 



cent are broad-linc and narrow-line AGN respectively. 
A small contamination (6 per cent) by stars is also ex- 
pected (Schmidt et al. 1998), but since they are ran- 
domly distributed over the sky their effect would be to 
dilute somewhat the meeisured correlation signal. From 
the work of Woods & Falilman (1997) we can deduce 
that for a stellar contamination of 6 per cent a reduc- 
tion in the observed correlation signal by ~13 per cent 
should be expected (correcting for this reduces our To 
values, derived in section 4 by only ^^5%) . In the ROSAT 
Lockman Hole Survey there is also a small fraction of 
galEixy groups. As these may be more strongly clustered 
than galaxies, and possibly AGN, they may increase 
marginally the overall signal. For this reason we have ex- 
cluded the nine extended sources, which we have found 
in our XMM observations, from the subsequent analysis. 

The differential X-ray source counts in the 0.5-2 
and 0.5-8 keV spectral bands are shown in Figure 1 and 
are compared with the best fit relations of Baldi et al. 
(2002; soft band) and Manners ct al. (2003; total band). 
In the 0.5-2 keV band there is good agreement between 
our results and the Baldi et al. (2002) double power-law 
best fit to the number counts. The Manners et al. (2003) 
best fit is derived for sources in the flux range /(0.5 — 
8keV) = 10"^^ - 8 X lO^^'^ergs"^ cm"^. Although our 
dN/dS is in good agreement with their results in the 
above flux range, at brighter fluxes the surface density 
of X-ray sources is lower than the extrapolated Manners 
et al. (2003) relation. This suggests that a double power- 
law is required to fit the 0.5-8 keV dN/dS over the fiux 
range 10~^^ — 10~^^ erg s~^ cm~^. We therefore adopt a 
double power-law of the form: 
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where A2 = Ai + (Si - B2) x fx, fx is. the flux 
at the break. We estimate /x(0.5 — SkeV) ^ 6 x 
10"" ergs"^ cm-^ Ai = -8.9 ± 2.2, Bi = -1.8 ± 0.2 
and B2 = —2.3 ± 0.1. Our best-fit double power-law 
relation is shown in Figure 1. 



3 TWO-POINT CORRELATION 
FUNCTION ANALYSIS 

The two-point angular correlation function, w{0), is de- 
fined as the joint probability of finding sources sepa- 
rated by an angle 6. For a random distribution of sources 

■w{6) = and therefore, the angular correlation function 
provides a measure of galaxy density excess over that ex- 
pected for a random distribution. In this paper wc use 
the estimator described by Efstathiou et al. (1991) 



(1) 



with the uncertainty in w{6) is estimated from the re- 
lation 



(Ty, = ^/{l+w{e))/NDR, 



(2) 



where No d is the number of data-data pairs in the in- 
terval [9 — AS, 9 + AS] and Ndr is the number of data- 
random pairs for a given separation. In the above rela- 
tion / is the normalization factor / = 2Nii/{ND — 1) 
with Nd and Nr being the total number of data and 
random points respectively. For each XMM pointing we 
produce 100 Monte Carlo random catalogues having the 
same number of points as the real data which also ac- 
count for the sensitivity variations across the surveyed 
area (see section 2). Furthermore since the fiux thresh- 
old for source detection depends on the off-axis angle 
from the center of each of the 'XM.'M- Newton pointing, 
the sensitivity maps are used to discard random points 
in less sensitive areas. This is accomplished by assigning 
a fiux to each random point using the differential source 
counts plotted in Figure 1. If that fiux is less than 5 
times the local rms noise at the position of the random 
point (assuming Poisson statistics for the background) 
this is excluded from the random data-set. Wc have ver- 
ified that our random simulations reproduce both the 
off-axis sensitivity of the detector as well as the individ- 
ual field logAf-logS. 

Using the methods described above we estimate 
w(9) in logarithmic intervals with 5\og9 ~ 0.05. For 
both samples we estimate w{9 < 150 ) ~ 0.11 ± 0.03 
corresponding to a statistically significant signal at the 
« 3.5(T confidence level (Poisson statistics). We now fit 
the measured correlation function assuming a power-law 
of the form w{9) = {9o/9)''-^, fixing 7 to 1.8. We use a 
standard minimization procedure: 



E 



wxMM(eO - {9o/9' 



17-1 



(3) 



with each point weighted by its error (cr*). Note, that the 
fitting is performed for angular separations in the range 
40-1000 arcsec. We also note that our results are insen- 
sitive to both the upper cutoff limit in 9 and the angular 
binning (for more than 10 bins) used to estimate w{9). 
Therefore, the best fit parameters for both the soft and 
the total band sub-samples are: 9o = 9.3 ± 1.9 and 9o = 
9.0 ± 1.7 arcsec's respectively. Note that the errors cor- 
respond to 1(7 (Ax^ = 1.00) uncertainties, which are es- 
timated using the vaxiation of Ax^ = X^(^o) — Xmin(^o) 
[Xmin is the absolute minimum value of the x^]- How- 
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Table 1. Angular correlation function analysis results. The columns are as: X-ray sub-sample, number of objects in the sub- 
sample, the corresponding angular correlation length, the reduced x^i the probabilities and the correlation signal between 
0-150 arcsec. The errors rcprosoiit la imrortaiutios. 
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0.5-8 kcV 


462 


10.8 ± 1.7 


1.50 


0.10 


0.114 ±0.037 


0.5-2 keV 


432 


10.4 ± 1.9 


1.10 


0.35 


0.105 ±0.035 



ever, these raw values should be corrected for two pos- 
sible bias presented below. 

3.1 Integral constraint 

When calculating the angular correlation function from 
a bounded region of solid angle Q, corresponding to the 
area of the observed field, the background projected lo- 
cal density of sources is Ns/il (where Ns is the number 
of objects brighter than a given flux limit). However, 
this is an overestimation of the true underlying mean 
surface density, because of the positive correlation be- 
tween galaxies at small separations, balanced by neg- 
ative values of w{d) at larger separations. This bias, 
known as the integral constraint, has the effect of re- 
ducing the amplitude of the correlation function by 

uJn = -^ j j w{e)dQ.idQ.'2. (4) 

Clearly, evaluating necessitates a priori knowledge 
of the angular correlation function. A tentative value 
of wn using a range of w^O) by varying within Icr our 
results is: ton ~ 0.01. 

Adding ojn to each bin of our raw w{9) the integral 
constraint has a small but not negligible effect on the 
estimated correlation lengths. Indeed, for the 0.5-2 keV 
band repeating the fittings using ujn ~ 0.01 we find 
6o — 10.7 ± 1.9 arcsec and 6o — 11.1 ± l.Tarcsec for the 
soft and the total band respectively. 

3.2 Amplification bias 

Another bias that may affect the measured angular cor- 
relation function of our X-ray sources is the amplifi- 
cation bias (e.g. Vikhlinin & Forman 1995). The orig- 
inal quantification of this effect can be traced back to 
Kaiser's (1984) work which showed that smoothing of 
the galaxy distribution using a Gaussian kernel with 
size similar or larger to the correlation length of the 
underlying galaxy distribution increases the correlation 
function of the resulting density peaks compared to that 
of the underlying galaxies. Furthermore, the larger the 
smoothing radius the higher the amplitude of the cor- 
relation function of the resulting density peaks. 

In the present analysis we are faced with a similar 
situation since the Point Spread Function (PSF) Full 
Width Half Maximum (FWHM) of the XMM-Newton 
detector is of the same order of magnitude (~ 6 arcsec) 
with the measured angular correlation length (6o ^11 
arcsec). X-ray Sources separated by less than ~ 6 arc- 
sec will be observed as a single object. This is in effect a 
smoothing process, similar to that of the Kaiser's study. 



with smoothing radius roughly equal to the XMM- 
Newton PSF size. Vikhlinin & Forman (1995) studied 
the clustering properties of X-ray sources detected on 
ROSAT archival data and found that their measured 
w{9) was severely affected by the amplification bias due 
to the large FWHM of the ROSAT PSF. In our case 
we expect significantly less problems since the XMM- 
Newton PSF size is smaller than that of the ROSAT 
detector. 

We quantify this effect using an approach that is 
similar to that of Vikhlinin & Forman (1995). These 
authors used the Soneira & Peebles (1978) algorithm 
to construct correlated point processes with a variety of 
in-build correlation amplitudes. Then using a smoothing 
window with the size of the ROSAT PSF they wore able 
to determine that their measured 'w{6) was artificially 
enhanced by a factor of ~ 2.85. 

Our method is based on the concept that the galaxy 
correlation function due to its power-law nature could 
be considered a fractal. Therefore one can shift the am- 
plitude oiw{6) at different scales keeping its slope fixed. 
This simplifies our study since we do not need to con- 
struct different correlated point process having a partic- 
ular correlation length. Any correlated point processes 
that is described by power-law with the required expo- 
nent can be scaled to have a specific correlation length 
one wishes. For our study we use the publically avail- 
able ACDM Hubble volume cluster distribution which 
has a well defined power-law correlation function with 
an exponent 7 ~ 1.8 (Frenk et al 2000). 

Lets assume that the angular correlation function 
of the model catalogue above has a correlation length 
6o,c while, the true (unaffected from the amplification 
bias) correlation length of the X-ray point sources is 
0o,x- We can translate the angular scale of the model 
correlation function to that of the XMM correlations 
by multiplying the former scale by the factor: 

f — ^o,x /^o,c • 

We can now simulate the effect of XMM PSF smooth- 
ing on the scaled model correlation function by merging 
all the model pairs with separations less than the PSF 
FWHM (i.e. ~ 6 arcsec for the XMM-Newton) and then 
fit the model angular correlation function to obtain the 
best fit angular correlation length-scale and compare it 
to that of the XMM point source data. 

However, since we do not know the value of 6o,x 
but it is rather the value that we seek to find from our 
analysis, we apply an iterative procedure by which we 
change the value of 9o,x, and thus of /, until the result- 
ing scaled model correlation function (i.e. after smooth- 
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Figure 2. The two-point angular correlation function for the 
soft (open points) and total (crosses) bands respectively. The 
dashed line represent the best-fit power law ■w{6) = (9o/Q)^'^ 
for the soft (0.5-2)keV band (see parameters in Table 1), 
while the continuous line represents the best fit ACDM 
(Oa = 0.7 and h = 0.7) model in the framework of Basilakos 
& Plionis (2001) biasing model. In the insert we present the 
measured correlation function in linear scaling. 



ing) has an angular correlation length equal to that of 
the raw XMM point source correlations. 

The previous analysis shows that our XMM-A'^euiton 
observations are only marginally affected by the ampli- 
fication bias. The true underlying angular correlation 
length of the X-ray population is overestimated by ~3- 
4 per cent. Therefore, we conclude that the corrected 
(free of the amplification bias) correlation length of the 
^yiyi-Newton soft X-ray sources is about 9o — 10.4 arc- 
sec. 

We validate the above procedure by recalculating 
the amplification bias for the case of ROSAT observa- 
tions. We adopt a correlation length (free from ampli- 
fication bias) of Oo.x ~ 4 arcsec (Vikhlinin & Forman 
1995) and an effective smoothing scale of ~ 20 arc- 
sec similar to the ROSAT PSF FWHM. Our procedure 
gives that the expected amplified correlation length of 
the ROSAT sources is ~ 10.5 arcsec a factor of ~ 2.7 
higher than the true value, in excellent agreement with 
the Vikhlinin & Forman (1995) analysis. We are there- 
fore confident that our method and results are robust. 



3.3 The final angular correlation length 

After taking into account the corrections described 
above, we present the corrected angular correlations in 
Figure 2 for the soft and total spectral bands. The best 
fit parameters for both sub-samples are presented in Ta- 
ble 1. Our results are higher than those of Vikhlinin et 
al. (1995) who derive the angular correlation function 



from ROSAT pointed observations which have a com- 
parable efi'ective flux limit with our XMM pointings. 
The above authors find an angular clustering length of 
9o ~ 10 arcsec, which reduces to 4 arcsec after correc- 
tion for the amplification bias. Comparison with the re- 
sults of Basilakos et al. (2004) shows that the hard band 
sources are more strongly clustered, at least on angular 
projection, {9a^'"^ ~ 22 ±9 for 7 = 1.8) compared to the 
soft band sources. 



4 THE SPATIAL CORRELATION LENGTH 
OF THE XMM SOFT SOURCES 

4.1 Inverting Limber's equation 

The spatial correlation function can be modeled as (de 
Zotti et al. 1990) 



e(r,z) = (r/ro)-^ X {1 + z)-^^+''> , 



(5) 



where e parametrizes the type of clustering evolution. 
If e = 7 — 3 (ie., e = —1.2 for 7 = 1.8), the clustering 
is constant in comoving coordinates (comoving cluster- 
ing) , which means that the amplitude of the correlation 
function remains fixed with redshift in comoving coor- 
dinates as the galaxy pair expands together with the 
Universal expansion. Alternatively, in the e = —3 model 
the clustering is constant in physical coordinates, while 
e = reflects the stable clustering model (eg. de Zotti 
et al. 1990). 

We can relate the amplitude 9o in two dimensions 
to the corresponding three dimensional one, Vo, using 
Limber's integral equation (cf. Peebles 1993). For ex- 
ample, in the case of a spatially fiat Universe, Limber 
equation can be written as 



w{9) 



(6) 



where (j){x) is the selection function (the probability that 
a source at a distance x is detected in the survey) and 
X is the proper distance related to the redshift through 



x(z) -- 

with 

E{z) 



dt 



Ho J, E{t) ' 



1/2 



(7) 



(8) 



(see Peebles 1993). The number of objects in the given 
survey with a solid angle ils and within the shell {z,z + 
dz) is: 



diV 
"dT 



(9) 



Therefore, combining the above system of equations, the 
expression for w{9) satisfies the form 



w{9) 



-'-L (tviit) ''W^^i, 



C(r,2)dii (10) 



Note that, the physical separation between two sources, 
separated by an angle 9 considering the small angle ap- 
proximation, is given by: 
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" (1 + ^) 

Using eq.(5) and eq.(lO) we obtain: 



(11) 



r2Ho 



1 dTV 



,iV dz J 



-3-6+7 



a;T-i(2) 



-dz 



where iJ^ = r(i)r(^)/r(i). 

In order to perform the inversion we still need to de- 
termine the source redshift distribution dN/dz. Since we 
have no unbiased redshift information for our sources we 
can resort to a measure of dN/dz using an estimate of 
their luminosity function. In flux-limited samples, there 
is a degradation of sampling as a function of distance 
from the observer (codified by the so called selection 
function). The latter also depends on the evolution of 
the source luminosity function but it is independent of 
the cosmological model, used in the derivation of the 
luminosity function. Thus for our X-ray sources the se- 
lection function can be written as: 



[x] = / ^{L^,z)dL 

'imi„(^) 



(13) 



where ^{Lx,z) is their redshift dependent luminosity 
function. In this work we used the soft band luminos- 
ity functions of Miyaji, Hasinger & Schmidt (2000) and 
of Boyle et al. (1993). We also use different models for 
the evolution of the soft X-ray sources: a pure lumi- 
nosity evolution (PLE) or the more realistic luminosity 
dependent density evolution (LDDE; Miyaji et al 2000). 
In Fig. 3 we present the expected redshift distributions 
of the soft X-ray sources for three different luminos- 
ity functions and evolution models. The LDDE model 
predicts a redshift distribution shifted to much larger 
redshifts with a median redshift of z ~ 1.2 (see also Ta- 
ble 2) comparing with both the Boyle et al (1993) and 
Miyaji et al (2000) luminosity functions with pure lu- 
minosity evolution. It is very interesting that the source 
redshift distribution of the ROSAT Lochman Deep field 
(Schmidt et al. 1998), albeit having a flux limit slightly 
lower than of our survey, traces quite well the LDDE 
predictions (see histogram in Fig. 3), a fact that sup- 
ports this luminosity function evolutionary model. To 
quantify this claim we have performed a l^st be- 
tween the observed and theoretical redshift distribu- 
tions and found that the probability of consistency is 
0.45 (xVdf = 0.97), < 10"^ (xVdf = 8.5) and 0.04 
(xVdf = 2.1) for the LDDE, PLE (Miyaji) and the 
PLE (Boyle) models, respectively. 



4.2 Results 

Using eq.(12), the LDDE luminosity evolution model 
and e — —1.2 we find within the concordance cosmo- 
logical model a soft band correlation length of ro — 
16.4 ± 1.3 Mpc. This value is comparable to that of 
Extremely Red Objects (EROs), luminous radio sources 
(Roche, Dunlop & Almaini 2003; Overzier et al. 2003; 
Rottgering et al. 2003) and hard X-ray sources (Basi- 
lakos et al. 2004) which are found to be in the range 
r-Q — 12 — 19 Mpc. It is interesting to mention that 
also radio sources which contain an AGN show strong 




Figure 3. The redshift selection function for three different 
luminosity function models: (a) Miyaji et al (2000) with PLE 
(dashed line), (b) Miyaji et al (2000) with LDDE (continuous 
line) and (c) Boyle et al. (1993) with PLE (dot-dashed line). 
The histogram corresponds to the distribution of the Schmidt 
et al (1998) X-ray sources of the ROSAT Lochman Deep 
Field. 



clustering (ro — 11 /i ^ Mpc), while the opposite is true 
for the case of radio sources showing no AGN activity 
(Magliocchetti et al. 2004). 

However, our previously derived ro value is signifi- 
cantly larger than those derived from optical AGN sur- 
veys: ro ~ 5.4-8.6 h-^ Mpc (Groom & Shanks 1996; La 
Franca et al. 1998; Akylas et al. 2000; Groom et al. 2002; 
Grazian et al. 2004) as well as from the recent X-ray se- 
lected sample of MuUis et al. (2004) who find ro ~ 7.4 

Mpc. We can push our inverted ro values to approx- 
imate closely the latter results only if we use the con- 
stant in physical coordinates clustering evolution model 
(e — —3), in which case we obtain ro — 7.5 ± 0.6 
Mpc, which is in excellent agreement with the MuUis 
et al. (2004) results. Note that earlier ROSAT soft X- 
ray clustering results of Garrera et al (1998) found a 
weaker clustering, with their upper limit of the linear 
clustering evolution model being marginally consistent 
with our e — —3 results. 

In Table 2, we list the values of the correlation 
length, ro, resulting from Limber's inversion for different 
luminosity function, evolution models as well as for dif- 
ferent cosmological models. We can attempt to disentan- 
gle the different sources of the apparent ro differences. 
Firstly, comparing the LDDE model between the Ein- 
stein de Sitter and the concordance (ilm = 1 — = 0.3) 
Gosmological models it becomes evident that the effect 
of moving from the former to the latter model increases 
by ~ 50% the value of ro (for both e cases). Note also 
that within the EdS cosmological model moving from 
the LDDE to the PLE luminosity evolution models de- 
creases the value of ro by ^ 20% for e = —1.2, while 
for the e — —3 case there is no significant difference 
between the two luminosity evolution models. There- 
fore, although we do not have the PLE luminosity model 
parameters for the concordance cosmological model we 
may expect similar changes as before, which implies 
that within this cosmological and luminosity evolution 
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Table 2. The soft X-ray sources correlation length (ro in Mpc) for different clustering models (e) and for the different 
luminosity functions and evolution models. The last column indicates the predicted median redshift, from the specific luminosity 
function used. The bold letters delineate the preferred cosmological model and the most updated luminosity function. 



LF 


Evol. Model 


(Cm, Ha) 


ro (e = -1.2) 


ro (e = -3) 


z 


Boyle 


No evol. 


(1,0) 


7.9±0.6 


5.4±0.4 


0.50 


Miyaji 


No evol. 


(1,0) 


6.5±0.5 


4.9±0.4 


0.37 


Boyle 


PLE 


(1,0) 


12. Oil. 


6.3 ±0.5 


0.92 


Miyaji 


PLE 


(1,0) 


8.8±0.7 


5.7 ±0.4 


0.58 


Miyaji 


LDDE 


(0.3,0.7) 


16.4 ± 1.3 


7.5 ±0.6 


1.19 


Miyaji 


LDDE 


(1,0) 


11.2 ±0.9 


5.0 ±0.4 


1.19 



(PLE) models we would obtain ro ~ 12.5 and 7 h ^ 
Mpc for the e = —1.2 and e — —3 models, respectively. 

Also note that the change of the luminosity func- 
tion model and thus of the redshift selection function, is 
always accompanied by a change of the median redshift 
of the corresponding redshift distribution. 

We can attempt to parametrize the different lumi- 
nosity model effects on the determination of ro by in- 
vestigating its dependence on the median redshift of the 
source distribution as well as on the cosmological model. 
To do so we use a parametrized, by the characteristic 
redshift z, analytical selection function, given by Baugh 
(1996): 



dN 2 
— ; — oc z exp 
dz 



3/2 



(14) 
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Figure 4. Comparison of the expected increase of ro as a 
function of the median redshift of the source redshift distri- 
bution for a concordance (continuous line) and an Einstein 
de Sitter model respectively, using the parametric selection 
function model of Baugh (1996). The results of Table 2 based 
on the Miyaji et al (2000) luminosity functions are plotted 
(filled circles for the concordance model and open for the 
Einstein de Sitter model). 



where z = y/2zc is the median redshift. Although this 
formula has been derived from the distribution of op- 
tical galaxies while the redshift distribution of X-ray 
sources maybe different we find that at least for the 
case of the Miyaji et al (2000) luminosity function pro- 
vides absolutely consistent results. For example, insert- 
ing eq. (14) in eq. (12) with z ~ 1.19 and e = —1.2 we 
find for the LDDE model ro ~ 16.3 ± 2.0 Mpc and 
ro — 11.0 ± 1.5 Mpc for the concordance and Ein- 
stein de Sitter models, respectively, which are in excel- 
lent agreement with the direct LDDE results (see Table 
2). Similar consistency is found also for the other mod- 
els presented in Table 2 (except for the models based 
on Boyle's luminosity function which is due to their sig- 
nificant contributions from very large redshifts) . In Fig. 
4, we show the dependence of the derived ro on the 
median redshift of the source distribution for the two 
different cosmologies and two clustering evolution mod- 
els. We also plot our direct results (using the different 
Miyaji et al. luminosity function models that provide 
different z) of Table 2. The excellent consistency is ev- 
ident which makes us confident of our results. Guided 
by Fig. 4 we can deduce that the PLE model of the 
Miyaji et al (2000) luminosity function would provide 
ro — 12 and 7 Mpc within the concordance cosmo- 
logical model for the e — —1.2 and e = — 3 clustering 
models respectively. 



5 THE XMM SOURCES COSMOLOGICAL 
BIAS 

Within the framework of linear biasing (cf. Kaiser 1984; 
Benson et al. 2000), the mass-tracer and dark-matter 
spatial correlations, at some redshift z, are related by: 

ar,z)^^DM{r,z)b\z) , (15) 

where b{z) is the bias evolution function. 

We can quantify the evolution of clustering with 
epoch presenting the spatial correlation function of the 
mass ^DM(r, 2) as the Fourier transform of the spatial 
power spectrum P{k): 



eDM(r,2) = (l+2)-^^+^^2^ 



■fc2p(fc)!M^dfc,(16) 
kr 



where k is the comoving wavenumber. Furthermore, 
the predicted spatial correlation function of the X-ray 
sources can be written as: 



■ 2 r,/, N sin(fcr) , , 
k P{k) — e — ^dfc 



f (r z)-^ 
where 

R{z) = {1 + zy^^'+'h^z) 



(17) 



(18) 



As for the power spectrum of our CDM models, 
we use P{k) ~ k^T'^{k) with scale-invariant (n — 1) 
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primeval inflationary fluctuations and T{k) the CDM 
transfer function. In particular, we use the transfer 
function parameterization as in Bardeen et al. (1986), 
with the corrections given approximately by Sugiyama 
(1995). Note that wc also use the non- linear corrections 
introduced by Peacock & Dodds (1994). 



5.1 Bias Evolution 

The concept of biasing between different classes of ex- 
tragalactic objects and the background matter distribu- 
tion was put forward by Kaiser (1984) and Bardeen et 
al. (1986) in order to explain the higher amplitude of 
the 2-poiiit correlation function of clusters of galaxies 
with respect to that of galaxies themselves. 

The deterministic and linear nature of biasing has 
been challenged (cf. Bagla 1998; Dekel & Lahav 1999) 
and indeed on small scales (r < 10/i~^Mpc) there are 
significant deviations from 6(r) — const. Despite this, 
the hnear biasing assumption is still a useful first order 
approximation which, due to its simplicity, it is used 
in most studies of large scale clustering (cf. Maglioc- 
chetti et al. 1999). In this paper however, we will work 
within the paradigm of linear and scalc-indcpcndont 
bias. Based on different assumptions a number of bias 
evolution models have been proposed (eg. Nusser & 
Davis 1994; Fry 1996; Mo & White 1996; Matarrcsc et 
al. 1997; Tegmark & Peebles 1998; Bagla 1998; Plionis & 
Basilakos 2001). However, here wo will discuss two that 
have been shown to describe relatively well the evolution 
even beyond ~ 1. 

• Merging Bias Model (hereafter B2): Mo & White 

(1996) have developed a model for the evolution of the 
the so-called correlation bias, using the Prcss-Schechter 
formalism. Utilizing a similar formalism, Matarrese et 
al. (1997) extended the Mo & White (1996) results 
to include the effects of different mass scales (see also 
Moscardini et al. 1998; Bagla 1998; Catelan et al. 1998; 
Magliocchetti et al. 2000). In this case we have that 



bB2{z) =0.41-1- 



(bo - 0.41) 



(19) 



with /? ~ 1.8. Note that D{z) is the linear growth rate 
of clustering (cf. Peebles 1993) * scaled to unity at the 
present time. 

• Bias from Linear Perturbation Theory (hereafter 
B3): Basilakos & Plionis (2001, 2003), using linear per- 
turbation theory and the Friedmann-Lemaitre solutions 
of the cosmological field equations have derived analyt- 
ically the functional form for the evolution of the linear 
bias factor, h, between the background matter and a 
mass-tracer fluctuation field. For the case of a spatially 
flat A cosmological model (fim + = 1), the bias evo- 
lution can be written as: 

6b3 {z) = AE{z) + CE{z)K{z) + 1 (20) 
with 

* D{z) = {1 + z)~^ for an Einstein-de Sitter Universe. 
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Figure 5. The Vciriation of Ax^ around the best bias fit (60) 
using different clustering behaviours (left panel for e = — 1.2 
and right panel for e = —3). Note that the solid and dashed 
lines represent the bias from the linear perturbation (B3) and 
the merging (B2) bias models respectively. 



K{z) 



POO 

Jl+z 



dy 



K{z) = {l + z) 



-1/2, 



13 7 Oa 

6'2'6' flmil + zy-i 



(21) 



(22) 



where F is the hyper-geometric function. Note that this 
approach gives a family of bias curves, due to the fact 
that it has two unknown parameters, (the integration 
constants A,C). Basilakos & Plionis (2001, 2003) com- 
pared the B3 bias evolution model with other models as 
well as with the HDF (Hubble Deep Field) biasing re- 
sults (Arnouts et al. 2002; Malioccietti 1999), and found 
a very good consistency. In this work, for simplicity, we 
fix the value of C being ~ 0.004, as was determined in 
Basilakos & Plionis (2003) from the 2dF galaxy corre- 
lation function. It is evident that the bias factor at the 
present time can be obtained from eq.(20) for = 



bB3{0)=A + CK{0) + l 

where we have used E{0) 
for Qa = 1 - fim = 0.7. 



(23) 

1. Note that K{0) ~ 9.567 



5.2 The bias at the present time 60- 

Based on the luminosity dependent density evolution 
(LDDE; Miyaji et al 2000) we quantify the bias factor 
at the present time 60, performing a standard min- 
imization procedure between the measured correlation 
function for the soft band (0.5-2)keV with that expected 
in our ACDM cosmological model, 



UIxMm(^') - Wmodel(^',6o) 



(24) 



where a* is the observed w{9) uncertainty. 

In Fig. 5 we present for the two bias and cluster- 
ing evolution models the variation of Ax^ ~ X^i^o) — 
Xmin(^o) around the best 60 fit, while in Table 3 wc list 
the results of the corresponding fits for all the consid- 
ered models. The resulting present time bias is between: 
60 ^ 1.05 - 1.90 and bo ~ 1.64 - 2.74 for the B2 and 
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Figure 6. The function R{z) = (1 + z)-(^+''h'^(z) as a 
function of redshift. The continuous hne {B3 bias) and the 
filled point (B2 bias) types represent the R{z) behaviour in 
the framework of a comoving clustering (e = —1.2), while 
the dot dashed line (B3 bias) and the crosses (B2 bias) are 
based on a clustering model which is constant in physical 
coordinates (e = —3). In the insert we present the difference, 
Ab{z) = bg^{z) — 6b2(-2)i between the {B3) and (B2) bias 
models as a function of z for e = —1.2 (continuous line) and 
e = —3 (dashed line) respectively. 



Table 3. Results of the predicted soft X-ray sources bias. 
The correspondence of the columns is as follows: bias and 
clustering evolution models, bo is the bias at the present time, 
the reduced a-nd the probabilities. Errors of the fitted 
parameters represent Icr uncertainties. 



Bias Model 


e 


bo 


xVdof 


P^2 


B2 


-1.2 


1.90 ±0.10 


0.97 


0.49 


B3 


-1.2 


2.74 ±0.15 


0.86 


0.60 


B2 


-3.0 


1.05 ±0.05 


0.97 


0.49 


B3 


-3.0 


1.64 ±0.05 


0.84 


0.63 



ing, e — —1.2 (continuous line for B3 and filled points 
for B2), the distribution of soft X-ray sources is strongly 
biased (1.90< bo< 2.74), as opposed to tiie less biased 
distribution (1.05^ bo<, 1.64) in the e — —3 (dashed line 
for B3 and crosses for B2) clustering model. This is to be 
expected, simply because the value e = — 3 removes the 
{1 + z) dependence from the R{z) functional form and 
thus, produces a lower corresponding correlation length 
(see Table 2), in contrast with the comoving (e = —1.2) 
clustering case. In other words, the higher or lower cor- 
relation length corresponds to a higher or lower bias 
at the present time respectively, being consistent with 
the hierarchical clustering scenario (cf. Magliocchetti et 
al. 1999). Note, that the above predictions are in good 
agreement with those derived by Treyer et al. (1998), 
Carrera et al. (1998), Barcons et al. (2000) and Boughn 
& Crittenden (2004), who have found feo ~ 1 — 2. 

Regarding the predictions of the two bias mod- 
els, we present in the insert of Fig. 6 the difference, 
Ab{z) — bB3{z) — 632(2), between the B3 bias model 
and the Mataresse et al. (1997) model B2 as a func- 
tion of redshift. The B2 bias evolves significantly more 
than the B3 model at relatively low redshifts {z < 1.0), 
which could be attributed to our assumption that the 
galaxy number density is conserved in time. It is evi- 
dent that merging processes, not taken into account in 
the B3 model, are probably important in the evolution 
of clustering. 

It is evident that the behaviour of the inverted X- 
ray source spatial correlation function is sensitive to the 
difi'erent values of e but there is also a strong dependence 
on the bias models that we have considered in our anal- 
ysis. 

We can attempt to select the most viable bias and 
the clustering evolution models by: 

(i) invoking the results of the local X-ray AGN clus- 
tering of Akylas et al (2000) and MuUis et al. (2004), 
who find To — 6.5 and 7.4 Mpc, respectively and 

(ii) noting that the local galaxy distribution, with a 
correlation length To — 5 Mpc, is unbiased with 
respect to the corresponding underline matter distribu- 
tion (eg. Lahav et al. 2002; Verde et al. 2002). 

These two facts leads us to a local bias between the X- 
ray selected AGN population and the underline matter 
distribution of: 



B3 bias models, respectively. Note that the theoretical 
(ACDM + B3 bias model) fit to the measured soft X- 
ray source angular correlation function is presented as 
the solid line in Fig. 2. 

In order to understand better the effects of AGN 
clustering, we present in Fig. 6 the quantity R{z) (see 
eq.18) as a function of redshift for the concordance cos- 
mological model and for different bias evolution models. 
It is quite obvious that the behaviour of the function 
R{z) characterizes the clustering evolution with epoch; 
in general AGN clustering is a monotonically increasing 
function of redshift for both B2 and B3 biasing mod- 
els. Figure 6, for example, clearly shows that the bias at 
high redshifts has different values in the different clus- 
tering models. In particular, for the comoving cluster- 



6(0) = (ro,m/ro,AGN)' 



(5/7)- 



1.35 



which is consistent only with the e = — 3 model of clus- 
tering evolution while it is in between the predictions of 
the two bias models used. 



6 CONCLUSIONS 

We have studied the angular clustering properties of the 
soft (0.5-e keV) X-ray point sources found in the XMM- 
Newton/2dF survey. We find that there is a strong 
(3ct) clustering signal. Indeed, if the two point angu- 
lar correlation function is modeled as a power law, 
w{d) = {60 /6)^'^, then after correcting for the integral 
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constraint and the amplification bias the best-fitting an- 
gular clustering length is 60 ~ 10.4 ± 1.9 arcsec. 

Inverting Limber's equation and using the preferred 
luminosity dependent density evolution model for the 
luminosity function gives ro ~ 16 and 7.5 h'^ Mpc, 
for the constant in comoving (e = —1.2) and in phys- 
ical (e = —3) coordinates clustering evolution models, 
respectively. In the former case, the values for the clus- 
tering length are comparable with those of Extremely 
Red Objects (EROs) and luminous radio sources, and 
are significantly higher than those found from previous 
ROS AT suTveys (e.g. VikhUnin et al. 1995, Akylas et al. 
2000, Carrora ct al. 1998; MuUis ct al. 2004) and opti- 
cal QSO surveys such as the 2QZ (Groom et al. 2002) 
and that of Grazian et al. (2004). However, we obtain 
a quite good agreement with the above surveys, only 
in the latter case of a clustering evolution model where 
the clustering length remains constant in physical coor- 
dinates (e — —3). 

Comparing the measured angular correlation func- 
tion for the soft band (0.5-2)keV X-ray sources with 
the theoretical predictions of the preferred ACDM cos- 
mological model {flm = 1 — ~ 0.3) and two bias 
evolution models, we find that the present bias values 
is in the range of 1.9<j bo<, 2.7 for the e = —1.2 model 
and 1.0^ 60^ 1-6 for the e = — 3 model. 
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